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Abstract. Spin currents may be generated by applying bias voltages V to 
the nanostructures even in the absence of spin-active ferromagnetic interfaces. 
Most theoretical proposals concentrate on a concrete spin-orbit interaction and 
on the disorder-averaged effect. It remains underappreciated that any spin-orbit 
interaction produces random spin currents with a typical amplitude ~ e 2 V/h not 
affected by disorder. This work addresses such mesoscopic fluctuations of spin 
currents for generic model of a nanostructure where several quantum connectors 
meet in a single node. The analysis is performed in the framework of recently 
developed quantum circuit theory of Gq corrections and reveals four distinct 
mechanisms of spin current fluctuations. The results are elaborated for simple 
models of tunnel and ballistic connectors. 
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1. Introduction 

In recent years considerable theoretical and experimental work is aimed at the 
controlled manipulation of electron spin in (nanoscale) solid state systems, a field 
commonly referred to as spintronics pQ . Spin current defined as a flow of electron spin 
is one of the central and most useful concepts of this field. [2] Much confusion can be 
easily brought by irresponsible and/or uncareful use of this concept. This confusion 
is related to two facts. Firstly, spin currents, unlike electric or particle currents, 
are not usually conserved under interesting circumstances. Secondly, it is frequently 
not evident that the spin currents due to electrons near Fermi surface really present 
a dominant contribution in spin balance: This is certainly not true for equilibrium 
and/or electrically isolated systems. 

Therefore, it seems prudent to start the article with a short description of 
the measurement that can directly access such spin current [3] and/or its temporal 
fluctuations [4] Let us consider a current injector of the resistance R that pushes spin- 
polarized electrons into a metallic lead. In principle, this polarization disappears far 
in the lead. However, spin relaxation is usually a slow process so the spin current 
is preserved at significant distance L s f from the injector. Let us set up a spin- 
active interface within this distance from the injector. A spin-active interface provides 
different resistance for electrons of two different spin projection on the polarization 
axis of the interface. Most common example is the interface between a ferromagnetic 
and non-ferromagnetic metal. In this case, the polarization axis coincides with the 
direction of magnetization of the ferromagnet. 

The voltage drop over the interface is thereby sensitive to spin of the electrons 
that traverse it. We would like the difference of the voltage drops at opposite directions 
of the polarization axis is directly proportional to the projection of the spin current 
on the polarization axis. This enables the measurement of spin current as a voltage 
measurement. To achieve this, one has to choose the interface resistance Ri n such 
that Ri n <C R but still by far exceeds the resistance of the metal piece between the 
injector and the interface. We shall stress the useful analogy between spin current and 
spin accumulation, on one side, and electric current and voltage on another side. The 
above condition on Ri n is precisely the same as the condition for the measuring the 
electric current from the injector through the voltage drop over the interface. 

The disadvantage of this measurement scheme is that at a given polarization of 
the interface only a single projection of the spin current can be measured. The point 
is that the spin current with the projection perpendicular to the polarization axis is 
usually immediately adsorbed at the interface. [5] An alternative scheme [5] where all 
three components of the spin current can be measured in the least obtrusive way is 
still difficult to realize for solid-state devices. 

Spin currents in the systems where spin-active interfaces dominate the resistance, 
are well-studied. [2] Usually, a simple circuit-theory approach based on the balance 
of spin and charge currents [6] suffices to describe everything. Here, we consider the 
situation where there are no spin-active interfaces in the nanostructure used as an 
injector of the spin current. Still there is spin-orbit interaction. 

It was noticed quite long ago [7] that in bulk homogeneous materials without 
inversion center electric currents produce spin currents. However, in bulk materials 
this effect is usually negligible in comparison with other agents of spin orientation, 
for instance, with the effect of magnetic field produced by the currents. [5] In ballistic 
nanostructures, the spin-orbit interaction in combination with geometric effects may 
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lead to spin currents. Consider, for instance, a single isotropic scattering center subject 
to electron flow in z direction. Intensity of electron scattering to +x and — x directions 
will depend on the spin projection on y axis. If we direct the scattered electrons to 
different terminals, we get spin currents in the terminals. Such device proposal has 
been actually elaborated. [9] At more abstract microscopic level, spin-orbit interaction 
always results in a spin-dependent scattering matrix (see, for instance, [10j). However, 
since the effect relies on geometry, it is drastically decreased by any scattering in the 
nanostructure. Indeed, the scattering makes electron distribution function isotropic 
so that they forget the geometry. In scattering matrix language, the spin currents 
are determined by random phase shifts of the scattering amplitudes. After averaging 
over random phase shifts, the average spin current is supposed to average out if the 
scattering is sufficiently chaotic. This has been confirmed by extensive numerical 
experiments. 

It remains underappreciated that even for completely chaotic scattering there 
must be significant mesoscopic fluctuations of spin currents. They are governed by 
the same quantum interference mechanism as the celebrated Universal Conductance 
Fluctuations [12] and must be of the same scale. Such fluctuations are clearly seen in 
numerical simulations at quantum level |13j . Recent numerical work |14j provides a 
detailed study of these fluctuations for a concrete device. 

In this work, we theoretically investigate the mesoscopic fluctuations of spin 
currents for multi-terminal nanostructures assuming sufficient isotropization of 
electron distribution function that allows to disregard all geometric effects. This 
means that after averaging over random phase shifts/impurity positions there is no net 
spin effect and the system is SU(2) invariant with respect to rotations in spin space. 
Under these conditions, spin-orbit interaction only manifest itself in spin relaxation 
and completely characterized by the rate of such relaxation. We denote the full spin 
current with spin index a — 1,2,3 coming to the terminal i as We study the 
correlators of the two, 

= Scflij 

that are SU(2) invariant. 

We stress that we study spin current fluctuations that persist at least at the 
scale of the contact size. This is very different from mesoscopic fluctuations of spin 
accumulation. Such fluctuations are known to form a random pattern correlated at 
the scale of mean free path [T5l [16] and are disregarded here. 

2. Model, parameters, expectations 

We model a nanostructure along the lines of quantum circuit theory [TTJ, [18] 
representing a complicated setup in terms of three kinds of the elements: terminals, 
nodes, and connectors. An arbitrary nanostructure with sufficiently chaotic scattering 
inside can be presented in this way with using sufficiently many nodes and connectors. 
The specific model under consideration consists of a single node attached to many 
reservoirs by arbitrary connectors (Fig. [1]). The node in circuit theory is very close to 
"chaotic cavity" used in Random Matrix Theory (RMT) of quantum transport. [19] 

(i) 

The connector i is characterized by a set of transmission eigenvalues T v , , p 
labelling quantum channels in the connector. Conductance of the connector reads 
£?(') = Gq J2 P ■ In the limits of applicability of the circuit theory, Gj) 3> Gq . 
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Figure 1. Generic model of nanostructure. We consider a multi-terminal setup 
(three terminals shown in the Figure). The connectors are characterized by the 
sets of transmission eigenvalues (Tp ). The node is characterized by Thouless 
energy related to the electron dwell time in the node. The spin-orbit interaction 
is represented by tso> spin relaxation time in the node. 

There is a single terminal per connector most generally characterized by an energy- 
dependent electron filling factor fi(E). The mesoscopic fluctuations of two-terminal 
conductance in this model have been considered in [20] by RMT methods. 

Spin current fluctuations can not be accessed by RMT since, as we will see, they 
disappear in the limit of pure ensembles. We need more parameters that are absent 
in RMT but present in more detailed and microscopic theories [12] that describe 
transitions between the ensembles. The size of the node is characterized by the 
mean level spacing Ss, or, more adequately, by the Thouless energy Eth = G^Ss /Gq, 
Gs == J2i G^j} . Mesoscopic fluctuations are supposed to correlate at energies ~ Eth- 
For the present purposes, the most important parameter is the spin-orbit time t$o 
that determines spin relaxation in the node. We assume no spin-orbit effects in 
the connectors. Therefore, t$o contains all information concerning the spin-orbit 
interaction and its strength is given by a dimensionless rjso = f^l^soEth- 

At moderate values of this parameter, r/so — 1, we expect the mesoscopic 
fluctuations of spin current to be very similar to the mesoscopic fluctuations of the 
charge current. Those are known to be universal corresponding to the fluctuating 
conductance of the order of Gq. At low voltages eV <C E t h this gives {(I 2 )) — GqV 2 
while at large voltages eV > E th {{I 2 )) ~ G 2 Q VE th /e [21]. The distinction is that 
the spin current is not conserved. Therefore, we expect non-zero fluctuations of the 
total spin current I a — '^2 i Ia^ to/from the nanostructure. Those must be of similar 
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magnitude as the fluctuations of the current to/from a reservoir. 

Naturally, we expect the spin current fluctuations to vanish in the limit rjso "C 1 
of weak spin-orbit interaction. The opposite limit is less obvious: One might think 
that strong spin-orbit interaction leads to fast spin relaxation in the node, the latter 
obviously suppresses any spin effects. 

The results derived and discussed in the article conform to these general 
expectations and give details of spin current fluctuations and their correlations for 
connectors of low and high transmission. 



3. Circuit-theory action 

Qualitatively, the model is most conveniently described by the action S that depends 
on the matrix Green's function of the node G and those of the reservoirs, (5W. All 
such matrices obey the constrain G 2 = 1. It is also known that the concrete matrix 
structure of the matrices does not have to be specified at this stage: most terms in 
the action do not depend on concrete physical realisation of such structure. [18] 
The action reads: 

S = Sj+S E +S so ; (1) 



2 ^ 

p 



' ( T (i) > 
In II -L-(G^G + GG {i) - 2) 



(2) 



S E = -^Tr[EG] (3) 

OS 

S so = — —Tv[Ga a Ga a ] (4) 

odsT so 

Here, <Sj gives a contribution of an individual connector that depends on the matrix 
Green's functions G , at its ends and the transmissions of the transport channels. 
Next term, Se, accounts for the finite size of the node and therefore for the de- 
correlation of the mesoscopic fluctuations ((A(E)A(E'))} upon increasing energy 
difference E — E' . The last term presents spin-orbit effects, a a being Pauli matrices 
in spin space. 

We make use of Keldysh technique, our action approach is similar to [55] The 
Green functions in the terminals at given energy E have the standard Keldysh 
structure 

* _ / 1 - 2fi 2fi \ , , 

Gi ~ \ 2-2/, -1 + 2/, ) (5) 

fi(E) being the energy-dependent filling factor in each reservoir, and are diagonal in 
spin index. The optimal point of the action gives the equilibrated Green function 
Gat in the node. It has the same structure as {5} with the balanced filling factor 
/n — J2i /*/ Si Gj). From now on, it is convenient for us not to distribute 
explicitly transport channels over the connectors. We achieve this by labelling all 
transport channels with i and ascribing a terminal with filling factor /, to each 
transport channel. In these notations, /at = Si^/i/Si^- 

To account for spin currents to/from each reservoir, we perform rotations with 
spin-dependent counting field \f m eacn terminal[23j, 

Giixf) = e^p(tx a a a T z /2)G l exp(~t X a <J a T z /2) (6) 
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In this case, the action will have a non-trivial saddle-point Gpj({xf} with the x~ 
dependent optimal value of the action. The expansion in Xa gives the momenta of 
low- frequency temporal fluctuations of spin currents. [5] 

However, such temporal fluctuations is not a problem of present interest. Firstly, 
we are interested in the time-averaged spin current only, so we only have to keep the 
first term in ^-expansion of each Green's function, 

GiiXi) ~ Gi+ 5G f , 6G* = 2x«a a ( ° _ 1 M (7) 

Secondly, this time-averaged current vanishes in the main order in G/Gq so we 
shall turn to the analysis of so-called GQ-corrections corresponding to mesoscopic 
fluctuations. The general setup for the evaluation of Gq corrections in the framework 
of circuit theory has been put forward in Ref.[I8] and, as it is accustom for fluctuations, 
involves quadratic expansion of the action in the vicinity of the saddle point. It goes 
as follows. We wish to access the correlation of the mesoscopic fluctuations at two 
different values of a parameter set: we call these values "black" and "white". For us, 
the parameter set in principle includes energy E, filling factors ft at this energy and 
counting fields xt ■ We find the saddle points G^,G^ corresponding for each set of 
parameters ( E w , f™, x™' a and E b , ff, Xi'"*)- We double the matrix structure so that 
the node Green function consists of four blocks: two diagonal ww,bb and two non- 
diagonal wb,bw. We give G a fluctuation g concentrated in the non-diagonal blocks. 
The constrain G 2 = I is satisfied up to the second order in g if we substitute 

G = G N +g- g 2 G N /2 

under the constrain 

gG N + G N g = 0. (8) 

If we introduce a compound index a composed of two check indices, the result of 
quadratic expansion can be written as 

6S = gtM ul g£»-r, a gtr. (9) 

To take proper care of the constrain l[5]). we consider matrix II a ^ defined trough the 
following relation: 

- (g - G gGo) -> 

\ (g bw - G b g bw G w ) , (10) 

the last equation makes white-black block separation explicit. We note that is 
a projector. It separates "bar" space on two subspaces where g either commutes or 
anti-commutes with Gat, and projects an arbitrary g onto anti-commuting subspace. 
The matrices M ab ,H ab define the fluctuation correction Sg q sought: 



5 0<l =lndet^M a5 J; (If) 

M-ah = n s5 M 55 n g5 + 5- ah - n g5 . (12) 

The quantity of our interest is the correlator of two spin currents per energy 
interval. Since a single current is given by a term linear in x, the correlator is given 
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by the term which is linear in both x w & n d X b i 

dS 

ij?mf k (E w ))) = Jim G ° (13) 

Eventually, this gives diffuson contribution to Sg q ■ The Cooperon contribution is 
obtained in similar manner with one of the Green's functions (G w or Gb) transposed. 
However, our particular model is invariant with respect to time reversal. Diffuson and 
Cooperon contributions to spin current fluctuations are therefore identical and it is 
enough to calculate only one of the two. 



4. Strategy and workout 

The calculations by the method described are straightforward but, admittedly, rather 
lengthy. The correct choice of strategy is vital for speedy evaluation. An obvious 
strategy is to follow [18] in detail: compute the fluctuation action by expanding near 
the x-dependent saddle point, get it as a function of x, go to the limit (|13p . However, 
such straightforward approach An alternative strategy is to expand the Green function 
near the "trivial" saddle-point solution x = 0. Since the x-dependent saddle point is 
close to trivial point, it will be sooner or later covered by such expansion. In fact, 
this approach is very close to the traditional Cooperon-Diffuson technique. [21] A 
disadvantage is that in this case one has to expand to at least third order in g and 
satisfying the constrain in this order is a headache. 

Here, we adopt a mixed strategy. We expand in the x-dependent point but keep 
X small. This amounts to expanding the matrix M from Eq. rT2]to the orders in x w ' b 
we need, 

M = M + M w + M b + M wb , (14) 

where 

Mw,b oc x w,b ; Mwb oc x w XW, 
Expanding the log in Eq. rT2]we observe that the term we need is given by 



S Gq = Tr 



M^Mo)- 1 - M w {M r 1 M b (M )- y (15) 



We also observe that elements of the matrix (Mo) -1 correspond to "common" 
diffuson propagators of the traditional mesoscopic fluctuation technique. Two terms 
in (115|) contain one and two diffuson propagators and in traditional theory describe 
respectively local and distant current correlations. 

An essential element of the strategy is to evaluate M in a convenient basis. We 
choose the basis in Keldysh-spin space in such a way that Gat(x) always remains 
diagonal. At x — this is achieved by the transformation 

P-iG N P=-T z ;P=(f N _ 1 ;). (16) 

separately in black and white sector. The basis rotates upon changing x w ' b i but this 
rotation is easy to take into account perturbatively. This basis is especially convenient 
to satisfy the anticommutation constrain |5]): It is satisfied by 

» = (s-0 + )' (17) 
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for both non-diagonal blocks. Owing to SU(2) symmetry in spin space, spin structure 
of the mesoscopic fluctuation in a given block(say, bw) can be most conveniently 
presented in terms of the singlet and three triplet components, g bw + g^"o- a . The 
matrix Mq does not contain elements that mix the components. For singlet-singlet 
elements we obtain 

(M )+- = ((Mo)-+)* = 

~ (j2 T ^j {l-i(E b -E w )/E th ); 

(M ) ++ = \j2 T ?f*f™- ( 18 ) 

i 

Here we introduce conveniently physical notations: filling factor drops over the 
transport channel i, f t = fi — /jy. For instance, the charge current per energy interval 
in this channel is given by (2ir)~ 1 Tifi. For triplet-triplet elements, r/so comes into 
play: 

(M )+_ = «M )-+)* = 



(19) 



Further calculations give M b ,w: M wb . We give expressions for this matrices skipping 
the trivial spin structure and spin index of Xi- We also skip ++ elements since latter 
do not contribute to the trace (|T5|l . 

M h + _ = M b _+ = 

^(l-Ti)xJ^ (20) 

and similar for M w . Here we simplify the answer by introducing counting field "drops" 
over each connector 

Xi =Xi~ Xn,Xn = TWTn , r (21) 

The above relations are similar to those for filling factors apart from the factor (l+rjso) 
that accounts for spin relaxation in the node. Without the factor, the structure 
guarantees the current conservation: If \i are the same in each terminal, all the 
"drops" x vanish. As to M wb , 



M- b --\T, T ^xT + ^ 



VsoXnXn', (22) 



Mf_ = M™\ = j£l?(l- T^JhXfi- (23) 

i 

Thereby we give the most general answer for the fluctuations of the spin current 
J2^( E ME b )))xtX h i = (2^)- 2 ( 

i,k 

A b A w Re(D t (D s -D t ))+ (24) 

A wb Re(D s -D t )+ (25) 

+FB wb (\D s \ 2 -\D t \ 2 ) (26) 

+TVS0X%XN(\D s \ 2 -\D t \ 2 )) (27) 
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where 



A btW = Y,T i {l-T i )$ w f b i ' w /Y,T i 



i i 



A wb = e Tf{i m b jfxffr/ E T * 



i 



D, 



l + i(E b -E w )/E Th ' 



1 



(28) 



l + i(E b - E w )/E Th + vso 



It is instructive to confront this with the corresponding expression for the fluctuations 
of the charge currents, 



where the definition of x hi constants A, B differs from (f2"Tj) by setting -qso = 0. 

The correlators are linear in both f w , f b as it should be since each (spin) current 
is linear in filling factors. 

5. Discussion of the general result 

We see that the spin current fluctuations are in fact contributed by three distinct 
mechanisms corresponding to three terms (fM|) . (|2~5|) and (|26[) . 

First mechanism is the only one that provides fluctuations for the system of tunnel 
connectors where all Ti <C 1. Inspecting its form, we find that the fluctuating spin 
currents in each transport channels are always proportional to the corresponding filling 
factor drops, and, apart from (1— Ti) factor, to the corresponding charge currents. So it 
looks like the fluctuating quantity is the polarization: The ratio of conductances of the 
channel, and not the spin/charge arriving to the transport channel. Moreover, the form 
of the term suggests that these polarization fluctuations are "globally" orchestrated: 
The polarization fluctuation is the same for all transport channels. Let us call the 
corresponding mechanism global polarization fluctuations(GPF) . Such "globality" is 
in agreement with the fact that in a common diffuson technique such two-diffuson 
diagramms describe mesoscopic correlations at large distances. We also see the same 
global correlation in corresponding contribution to the charge current fluctuations: 
in this case, it looks like the conductances of all transport channels fluctuate all 
together by the same value. Since fast spin relaxation in the node destroys the "global" 
coordination, the mechanism ceases to work in the limit of large rjso- 

The polarization fluctuations do not have to be correlated at global level: One 
can easy imagine that they arise locally and do not correlate in different transport 



Y,{j™{E w )j b {E b )))xtx\ = {^)- 2 { 



i,k 




(29) 
(30) 
(31) 



A wb Re (D s + 3A) + 
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channels. This is how the second mechanism (|2"5j) works. We see that also in this 
case the current fluctuation in each channel is proportional to the corresponding filling 
factor drops. This proves that the fluctuation is a change of the polarization. However, 
we also see that in contrast to (pM|) . (1251) is contributed by the products of counting 
field drops at the same channel only, XiXi- Therefore, they do not correlate in different 
channels. Let us call it local polarization fluctuations(LPF) . Due to the locality, the 
spin current fluctuations survive in the limit of large rjso- Similar mechanism (Eq. 
[30]) provides uncorrelated fluctuations of channel conductance for charge currents. As 
to Ti dependence of the mechanism, the origin of (1 — Ti) factor is clear: The channel 
with Ti = 1 has the maximum possible conductance for both spin directions, so the 
fluctuations of conductance/polarization must be suppressed. The origin of Tf factor 
is not clear for the author. However, the factor is here and owing to this the LPF 
mechanism does not work in the limit of tunnel junctions. 

A note on locality is required at this place. Locality of the polarization 
fluctuations (seen as XiXi terms in the general answer) does not generally imply the 
locality of the current fluctuations (that would be seen as xTXi terms). The reason 
for this is evident: (Incomplete) spin current conservation in the node. Let us give a 
polarization fluctuation to a connector. It produces not only a spin current that goes 
to the terminal. At the same time, it sends a spin current (equal in magnitude and 
opposite in spin direction) to the node. The node re-distributes this current over all 
outgoing channels. Therefore, a local polarization fluctuation creates the spin current 
in all transport channels. If spin relaxation in the node is sufficiently strong, the 
re-distribution ceases (\n — * 0) and spin currents are really local. 

The third mechanism also does not work for tunnel junctions. However, it is 
one of the two that work for purely ballistic connectors (T = 1) that cannot exhibit 
conductance/polarization fluctuations. Inspecting (|26p . we note that the spin currents 
in each connector are no more proportional to local filling factor drops. Rather, their 
intensity is set by a single global factor F that depends on filling factor drops on all 
channels and defines. However, the fluctuations are local as for LPF mechanism, with 
the same reservation. The picture behind it is as follows. Let us note that the node 
in our model works like a distributor: It receives electrons coming through transport 
channels and sends then back. If it does so with a slight random preference of the 
channels with respect to spin, we reproduce the structure of (|25]) . Let us call the 
mechanism local distribution fluctuations(DF) . Similar to LPF, it survives the limit of 
large -q S o- 

The fourth mechanism l|27p does not manifest itself in charge current fluctuations 
since it is due to incompleteness of spin conservation in the node. As in the case of 
LDF, the intensity of current fluctuations in all channels is set by the global parameter 
F. In distinction from LDF, the fluctuations are globally orchestrated. It looks like our 
node-distributor has some supply of spin of random direction and sends it uniformly 
to all channels. Let us call it global distribution fluctuations(GDF) . The contribution 
of this mechanism vanishes upon increasing rjso- 

6. Tunnel Connectors 

We turn to more concrete and quantitative examples. Let us assume that all 
connectors are tunnel junctions, so that all Tj < 1. As mentioned, the fluctuations 
in this case are due to GPF mechanism. It is convenient to turn back from 
transport channels to terminals, so i now labels the terminals with conductance 
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Figure 2. Tunnel connectors a. Variations of charge (dashed line) and spin 
currents versus the strength of spin-orbit interaction nso- The variations are 
normalized by the value of charge current fluctuations at rjso = 0- Note factor of 
8 in the scale for spin current fluctuations, b. Ratio of variations changes slowly 
in a wide interval of rjso- 



Gi = Gq ^2 p Tp . Let us restrict ourselves to voltage differences <S eExh ■ Under 
these conditions, we may disregard energy differences in the diffuson propagators D s t . 

(i) 

Integrating the general answer over energies and taking the limit Tp — > 0, we obtain 
the following relation for the correlator of the spin currents: 

- 6 a /4 (p)' (_J 7TZ7, — ^ 

J \Gt,J \l + Vso {1 + VsorJ 

Here Ii is the average current to the terminal i determined by the conductances of 
the connectors and the voltages applied in accordance with Ohm's law. We can now 
explicitly see that the spin current fluctuation in each connector is proportional to 
the current, so that the fluctuating quantity is the polarization. The full spin current 
to/from the system does not fluctuate, although it is generally allowed. This is specific 
for tunnel connectors. 

It makes sense to compare the spin current fluctuations with fluctuations of charge 
currents. Those are given by 

«^»=Ki) 2 ( 1+3 (dsF) (33) 

As known, [121 119j this fluctuation reduces by a factor of 4 upon increasing rjso ■ I n 
RMT language, this is the transition from orthogonal to symplectic ensemble. We 
plot in Fig. 2 the strength of spin and charge fluctuations and its ratio versus rjso- 
Spin fluctuations reach maximum at rjso ~ 0.75. The ratio of the fluctuations reaches 
maximum at 7750 ~ 1-5 and falls off very slowly upon increasing 770- The polarization 
of the fluctuating current in maximum is about 50 %. 

7. Ballistic Connectors 



Let us turn to the analysis of purely ballistic connectors (Tj = 1). As in previous 
Section, we arrange transport channels into connectors and label them with i. The 
fluctuations are contributed by LDF and GDF mechanisms. Again we assume 
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Vi <C Exh- We integrate the non-equilibrium factor F over energies to obtain (apart 
from e 2 factor) a strikingly simple expression: 

F = Y,hVi/G* (34) 

i 

where Vi is the voltage drop at the corresponding connector. The non-equilibrium 
factor is thus proportional to the total energy dissipation in the system. 

The factor determines the strength of the spin current fluctuations. From the 
general answer (|2"r?)) we obtain the contribution of LDF mechanism: 

Wf» - ^ = ^-F (gA 3 - r !ff s )2 ) ^( 2 + ^) (35) 
J e 2 G s V Gz[l + r]sorJ (1 + Vso) 2 

We can explicitly see that the fluctuations are local ( with the reservation made above). 
It is instructive to give the expression for the charge current fluctuations: 

«/,/,» - ** - F (a*, - (l + ITT ^ p ) (36) 

The contribution of the GPF mechanism reads: 

//rarfn i ~ GjGj rj 2 so (2 + T] SO ) 

■> 8Gs(l + vsor [l + Vsor 

In distinction from the tunnel junction system, there is a fluctuation of the full 
current to/from the system, Ip = If . It simply reads 

win = ^ s f + ^ (38) 

J e z (l + r/sor 
where we sum up the contributions of both mechanisms. We plot in Fig. 3 the 
fluctuations of electric current, full and transport spin current for a simple system of 
two connectors with equal conductances. 



8. Conclusions 



To conclude, we have studied the mesoscopic fluctuations of spin currents for generic 
model of a single-node, multi-terminal nanostructure. Spin-orbit interaction is 
characterized by a single parameter rj$o which is the ratio of the dwell time h/ ' Exh in 
the node to spin relaxation rate tso- At moderate values of this parameter, the scale 
of spin current fluctuations is the same as that of charge current fluctuations. 

We have found that the spin current fluctuations are contributed by four distinct 
mechanisms: Fluctuations of polarization or distribution, both can be either local or 
global. Local mechanisms survive in the limit of strong spin-orbit interaction. If all 
connectors are of tunnel nature, only GPF works. For purely ballistic connectors, 
LDF and GDF contribute. We present simple formulas valid for tunnel and ballistic 
connectors and some quantitative results. 
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Figure 3. Ballistic connectors. We concentrate on a simple circuit with identical 
connectors (on the left). We plot variances of charge current (a), spin current 
in a connector(b), and full spin current to/from the node. The variances are 
normalized on the charge current variance at r)$o- 



References 

[1] Zutic I, Fabian J and Das Sarma S 2004 Rev. Mod. Phys 76 323 

[2] Brataas A, Bauer GEW, and Kelly PJ 2006 Phys. Rep. 427 157 

[3] Jcdema FJ, Filip AT, and van Wees BJ 2001 Nature 410 345; Schmidt G, Ferrand D, Molenkamp 

LW, Filip AT, and van Wees BJ 2000 Phys. Rev. B 62 R4790 

[4] Lorenzo AD and Nazarov YV 2005 Phys. Rev. Lett. 94 210601 

[5] Di Lorenzo A and Nazarov YV 2004 Phys. Rev. Lett. 93 046601 

[6] Brataas A, Nazarov YV, Bauer GEW 2000 Phys. Rev. Lett. 84 2481 

[7] Levitov LS, Nazarov YV, Eliashberg GM 1985 Zh. Eksp. Teor. Fiz. 88 229 

[8] Nazarov YV 1985 Fiz. Tv. Tela 27 1753 

[9] Kisclcv AA and Kim KW 2001 Appl. Phys. Lett. 78 775 

[10] Pareek TP 2004 Phys. Rev. Lett. 92 076601 

[11] Yamamoto M, Ohtsuki T, and Kramer B 2005 Phys. Rev. B 72 115321 

[12] Altshuler BL 1985 JETP Lett. 41 648 

[13] Kramer B, Dittmer K, Debald S, Ohe J, Cavaliere F, and Sassetti M 2004 Mat. Sci. -Poland 22 
445 

[14] Ren W, Qiao ZH , Wang J , Sun QF , and Guo H 2006 Phys. Rev. Lett. 97 066603 

[15] Nazarov YV 1990 JETP Lett. 51 426 

[16] Zyuzin AY and Serota RA 1992 Phys. Rev. B 45 12094 

[17] Nazarov YV 1999 Superlattices and Micro structures 25 1221 

[18] Campagnano G and Nazarov YV 2006 Phys. Rev. B 74 125307 

[19] Beenakker CWJ 1997 Rev. Mod. Phys. 69 731 

[20] Brouwer PW and Beenakker CWJ 1996 J. Math. Phys. 37 4904 

[21] Larkin AI and Khmelnitskii DE 1986 Zh. Eksp. Theor. Phys. 91 1815 

[22] Feigel'man MV, Larkin AI, and Skvortsov MA 2000 Phys. Rev. B 61 12361 

[23] Nazarov YV and Bagrets DA 2002 Phys. Rev. Lett. 88 196801 



